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ABSTRACT 

We compare cosmic microwave background lensing convergence maps derived from South Pole Telescope 
(SPT) data with galaxy survey data from the Blanco Cosmology Survey, the Wide-field Infrared Survey Ex- 
plorer, and a new large Spitzer/IRAC field designed to overlap with the SPT survey. Using optical and infrared 
catalogs covering between 17 and 68 deg 2 of sky, we detect correlation between the SPT convergence maps and 
each of the galaxy density maps at > 4a, with zero correlation robustly ruled out in all cases. The amplitude 
and shape of the cross-power spectra are in good agreement with theoretical expectations and the measured 
galaxy bias is consistent with previous work. The detections reported here utilize a small fraction of the full 
2500 deg 2 SPT survey data and serve as both a proof of principle of the technique and an illustration of the 
potential of this emerging cosmological probe. 
Subject headings: galaxies: structure — cosmic background radiation 
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1. INTRODUCTION 

Gravitational lensing of the primordial anisotropies of the 
cosmic microwave background (CMB) imprints information 
about the density fluctuations between z ~ 1100 and the 
present day onto the observed CMB fluctuations. This infor- 
mation can be extracted by measuring the induced correlation 
between initially independent spatial modes of the CMB and 
used to construct a map of the lensing convergence field. This 
field is closely related to the projected gravitational potential. 

The statistics of CMB lensing ma ps provide powerful con- 
straints on cosmologi cal parameters (Lesgourgues et al. 2006; 
de Putter et al. 12009] ). Such maps can also be combined with 
other tracers of Targe- scale structure to test cosmological mod- 
els and constrain properties of the tracer population. One such 
analysis involves cross-correlating lensing maps with galaxy 
catalogs. The correlation measures the average lensing sig- 
nal from the dark matter halos that host the galaxies and can 
be used to determine the halo bias. These measurements test 
models of the time evolution of co smic densi ty fluctuations 
and of primordial non-Gaussianity (Dalai et al. 2008 ; Jeong 
|et al.|2009D . 

CMB lensing is a young field. The first detections — using 
lensing-g alaxy cross-co rrelation s — were reported relatively 
recently ( [Smith et a l. 2007; |Hirata et al.|20 08). The first de- 
tections of lensing in CMB data alone have come from high- 
resolution CMB experiments, bo th through non-Gaussianity 
( |Das et al.|201 lb| [van Engelen et al. 2 0T2]) and through smear- 
ing of the acoustic peak s (Reichardt et al. 2009 1 |Das et al. 



201 la[ |Keisler et al. 2011 ). Within the n ext year, further ad 
vances are expected: the Planck satellite (The Planck Collab 



oration 2006) will create all- sky convergence maps ( |Hanson 



et al.||2blT|), and the now-completed 2500 deg 2 South Pole 



Telescope (SPT) survey (Carlstr om et al.|20TT) will produce 
complementary lensing maps that have significantly higher 
signal-to-noise per mode than the Planck maps. 

Here we present the results of the cross-correlation of sev- 
eral galaxy populations with convergence maps from the SPT. 
In contrast to previous detections of this cross-correlation that 
utilized a larg e fraction of the sky with low signal-to-noise 
lensing maps ( |Smith et aLl|2007| |Hirata et aL|2008| ), we use 
convergence maps that have signal-to-noise greater than 1 on 
degree scales but significantly less sky area; the results pre- 
sented are from two ~ 100 deg 2 SPT fields. A significant 
cross-correlation is detected between the convergence maps 
and maps of galaxy density constructed from optical and in- 
frared (IR) catalogs in each of the two fields. 

The paper is structured as follows. We first discuss the un- 
derlying theory and provide a brief overview of the process of 
making convergence maps from CMB data. We next describe 
the catalogs, real and simulated, that we cross-correlate with 
the lensing maps. We conclude with a discussion of the re- 
sults and the potential of upcoming large CMB lensing data 
sets and large-area galaxy surveys. 

2. THEORY 

The CMB lensing convergence in a direction n on the sky 
is given in terms of the matter fluctuations by a line-of- sight 
integration, 

K(n) = Jd X W K (x)6(xn,z(x)), (D 

where S(r,z) is the fractional matter over-density at comov- 
ing position r and redshift z, and the distance kernel is 



( |Cooray & Hu|2000[ [Song et al.|2003) 
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Here, Q m is the matter density relative to the critical density 
evaluated today, Ho is the Hubble parameter today, a(x) is 
the cosmological scale factor, Xcmb — 14 Gpc is the comov- 
ing distance to the CMB recombination surface, and we have 
assumed a spatially flat Universe. 

Under the assumption that observed galaxies are biased 
tracers of mass fluctuations, the observed fractional galaxy 
over-density in a direction n is given by 
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where the distance kernel is 
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Here, dN(z)/dz is the distribution of galaxies in redshift and 
b(x) is the bias of the galaxies relative to the dark matter den- 
sity, assumed here to be independent of scale. 

The cross power between the convergence and the galaxy 
over-density at a multipole L is given in the Limber approxi- 



mation ( |Limber| 1 953} |Kaiser| 1 992] ) by 



tn> 
;K 



Cl 9 = j 



dz^^W K (x)W^ X )P(k 



(5) 



where the matter power spectrum today, P(fc, 0) is given by 
(<J*(k)<$(k')> = (27r) 3 (5 Dirac (k - k')P(fc, 0) under the assump- 
tion of independent Fourier modes k and k'. The amount 
of cross-correlation is thus determined by the overlap be- 
tween the two kernels, weighted by the matter power spec- 
trum. The factors which determine this overlap, namely 
1 lfb.«l (z)(P(k, z)dx/dz) 1 / 2 i are plotted for the CMB lens- 
ing convergence and several galaxy populations in the bottom 
panel of Figure [T] 

3. CONVERGENCE MAPS FROM CMB LENSING 

CMB maps of two fields, one centered at (RA,DEC) = 
(23h30m,-55d) and the other at (5h30m,-55d), were con- 
structed using 150 GHz data from the SPT survey. These 
fields together encompass 1 85 deg 2 and are the deepest SPT 
fields to date as they were observed with roughly twice the 
time per unit area as the rest of the SPT survey. I n addition 
to the 2008 data used in van Engele n et al. ( |2012| ) (hereafter 
VI 2), the maps used in this analysis include data from the 
2010-2011 observing seasons. As such, the resulting lensing 
maps are ~ 10% lower in noise than the typical maps to be 
expected from the rest of the survey. 

Conve rgence maps wer e constructed using a quadratic es- 
timator flHu & Okamoto||2002] ), as outlined in V12. Point 
sources with signal-to-noise greater than 6 were masked, 
while clusters detected via the thermal Sunyaev-Zeldovich 
(tSZ) effect were masked if they had a signal-to-noise greater 
than 6 in the 2008 data alone. This cross-correlation study 
is less sensitive to foreground contamination and offsets in 
power than the VI 2 analysis; therefore, the temperature maps 
were filtered to include modes from I = 1200-4000, a larger 
range than the I = 1200-3000 used in V12. 



The effect of foreground contamination is expected to be 
small for this analysis. Residual foreground emission in the 
CMB map that is correlated with the galaxy distribution is 
likely to lead to negative cross-correlations on large scales: 
excess small-scale power fluctuations in the temperature map 
are identified by the quadratic estimator as large-scale under- 
dense regions which 'squeeze' the sky. This produces an anti- 
correlation of true large scale structure with residual fore- 
ground contamination in the lensing reconstruction. Con- 
versely, dust in our Galaxy could also be interpreted by the 
lensing estimator as large-scale under-dense regions while 
also suppressing the galaxy counts, leading to a net positive 
correlation. We discuss tests for foreground contamination in 
Section 6. 

The measured CMB lensing signal was calibrated with sim- 
ulations. As outlined in VI 2, simulated lensed CMB skies 
including contributions from faint IR sources and tSZ from 
low-mass galaxy clusters (both modeled as Gaussian random 
fields) were added to noise generated from SPT timestreams. 
The simulations were filtered and analyzed in the same way 
as the real data. We normalize the reconstructed convergence 
maps as a function of L by comparing the cross-spectra of 
the input and output maps to the power spectra of the input 
maps. The simulated lensing maps are also used to estimate 
the uncertainty in the cross-correlation. 

4. OPTICAL AND IR-SELECTED CATALOGS 

We use optical- and IR-selected catalogs to create maps of 
the fractional galaxy over-density, (N — N)/N, where N is 
the galaxy counts in a 1 ; cell. A bright magnitude thresh- 
old is applied to each catalog to reduce contamination from 
Galactic objects, while a faint magnitude limit is set to ensure 
uniformity in source detection. We exclude regions around 
bright stars found in the 2MASS source catalog (Skrutskie 
|et al.|2006p . 

Optical catalogs were produced using data from the Blanco 
Cosmology Survey (BCS, (Des ai et aLl [2012)). The survey 
consists of two fields, a 5h and a 23h field, both contained 
within the SPT footprint. For this analysis, we extracted 
sources from images reduced using the pipeline described in 
|Restetal.| ( [200 5) and Mikna itis et"aL] ( |2QQ7b. Source detection 
was p erformed using SExtractor v2.8.6 ( jBertin & Arnouts 
1996), and the photometry was calibrated using Stellar Lo 



cus Regression (High et al. 2009 ). We restrict the analysis to 
areas that are complete to at least %ab = 22.5 and construct 
catalogs from all sources 19.5 < %ab < 22.5. 

In the 23h field, we also use data from the SPT Spitzer Deep 
Field (Ashby, private communication), an ongoing survey of 
100 d eg 2 at 3.6 and 4.5 /mi . Data reduction followed the meth- 
ods of As hby et al.| ( |2009] ), except, owing to the new survey's 
large size, photometry from individual 2 deg 2 sub-fields was 
merged in catalog space instead of extracted from a single 
monolithic mosaic. Here we analyze 32 deg 2 of existing cov- 
erage and include all sources with 4.5/im magnitude between 
15 and 17 Vega. 

As a second catalog in the 5h field we use the sources from 
the Wide-field Infrared Survey Explorer (WISE) (Wright et al. 



2010} all-sky data release. We restrict this analysis to areas 
Wl 



with Wl band exposure times between 231-616s, exclude re- 
gions with potentially corrupted objects as identified by the 
WISE pipeline, and include all sources with Wl magnitude 
between 15 and 17 Vega. Although the WISE selection func- 
tion is complicated at these magnitudes (owing to source con- 
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Fig. 1. — Redshift distribution for the galaxy samples are shown in the 
top panel. Mock catalogs (Section 5) are used to estimate the distribution 
for the optically-selected galaxies, while IR-selected galaxy distributions are 
based on photometric redshift estimates from Brodwin et al. (2006). The 
mock catalogs extend to z = 1.3, a sufficient redshift limit for the optical 
catalogs under consideration. The vertical dashed line at z = 1.5 indicates 
the maximum redshift of the Brodwin et al. (2006) spectroscopic training 
sample. The quantities shown in the bottom panel, which includes the CMB 
lensing power (green lines), are the curves that one would multiply to obtain 
dC 1 ^ 9 /dz, the integrand of equation 5, at L = 300 (solid curves) and L = 
1000 (dotted curves), assuming galaxy bias b(z) = 1. The dotted curves 
have been multiplied by 3 for clarity. 

fusionjand spatial variability in the imaging depth), we find it 
necessary to use these faint sources to obtain sufficient num- 
bers of high-redshift galaxies. Despite its limitations, the 
WISE sample is a particularly interesting test case, as WISE 
is an all-sky survey, and this analysis will soon be possible 
over the entire SPT survey area. 

The expected redshift distributions for the galaxy catalogs 
are shown in Figure [T] Summary statistics are presented in 
Table 1, and galaxy density maps are shown in Figure [2] 

5. MOCK CATALOGS 

Mock catalogs were used to estimate both the expected 
level of cross-correlation between CMB convergence maps 
and galaxy catalogs and the potential contamination from 
tSZ. The simulated galaxy sample is drawn from a 220 
deg 2 lightcone populated with galaxies in the redshift range 
0—1.3. The underlying dark matter distribution is based 
on a cosmological simulation of 1 /i _1 Gpc; this simula- 
tion is a single 'Carmen' simulation from the Large Suite 
of Dark Matter Simulations project (LASDAMAS) (see 
http://lss.phy.vanderbilt.edu/lasdamas/). The Adding Density 
Determined GAlaxies to Lightcone Simulations (ADDGALS) 
algorithm is run to assign galaxies to the dark matter particles 
in a way that reproduces the known luminosities, colors, and 
clustering properties of galaxies. The simulated catalogs have 
previo usly been used for tests of cluster finding (Koester et al. 
2007| ), photometric redshifts (Gerdes et al. 2010), and spec- 
troscopic followup (Cunha et al. 2011). A full description of 
the algorithm and the simulated sky catalog will be presented 
in Wechsler et al. and Busha et al. (in preparation). The 
simulated CMB lensing map based on the same underlying 
dark matter distribution was produced using multiple-plane 

1 http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/sec6_5.html 
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Fig. 2. — The top four panels show fractional galaxy over-density maps smoothed to 6' resolution. The bottom two panels show the lensing maps smoothed to 
1 ° resolution. At this scale, the signal-to-noise in the lensing maps is slightly greater than 1 . While the lensing maps shown here have been smoothed to highlight 
real mass fluctuations, it is clear from the cross-spectra there is substantial information at smaller scales as well. The color scale is ± 80% of the maximum 
deviation from the mean in each map. 

in the IR-selected catalogs. Also, these mock catalogs do 
not include AGN (which could enhance the cross-correlation), 
stars (which would not be correlated with a CMB lensing 
map), or the effects of source confusion (which would depress 
the cross-correlation). 

The simulated galaxy catalogs were constructed using the 
same selection criteria used for the real optical catalog. Maps 
of the fractional galaxy over-density were created and cross- 
correlated with the associated lensing convergence map. 



ray tracing methods similar to those of Hilbert et al. (2009), 
with modifications to produce maps with the proper lensing 
kernel at the redshift o f the CMB. Thermal SZ was included 
in the simulation as in Biesiadzinski et al. (2012): tSZ profiles 
were genera ted at the location of the halos using the |Arnaud| 
et al. ( 2010| ) model and projected along the line-of-sight to 
create a simulated tSZ map. 

Owing to the finite size of the simulation box, the mock cat- 
alogs and the CMB lensing map only extend to z ~ 1.3. For 
optical catalogs, this is not a serious limitation but there could 
be substantial cross-correlation coming from higher redshifts 



6. RESULTS AND DISCUSSION 



TABLE 1 
Field parameters and correlation statistics 



Field 



Area (deg ) Density ( 



sources 
deg 2 



) A(C L xlO~ 7 ) 



n 



X 2 (Best fit) Ax 2 (0) Bias 



WISE (5h) 


68.1 


BCS (5h) 


27.0 


BCS (23h) 


16.9 


Spitzer (23h) 


29.8 



6.9xl0 3 

2.5 xlO 4 

2.35 xlO 4 

1.4xl0 4 



0.19 ±0.05 
0.27 ±0.06 
0.24 ±0.07 
0.33 ±0.07 



-1.2 ±0.3 
-1.8 ±0.3 
-1.7 ±0.3 
-1.6 ±0.2 



11.3 
9.6 

13.7 



19.6 


0.9 ±0.2 


23.5 


1.2 ±0.3 


17.5 


1.1 ±0.3 


28.9 


1.7 ±0.3 



Galaxy catalog properties and results of power-law fits to lensing-galaxy cross spectra. For each catalog, we report best-fit amplitude and power law index, x 2 
of the best fit, the difference in x 2 from best-fit to a model with zero cross-correlation, and the bias. Note: we report the weighted area of the galaxy density 
maps multiplied by the lensing apodization mask. 
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Fig. 3. — Cross-spectra of galaxy number density maps with SPT lensing 
convergence maps. Spectra are calculated in 14 bins from L = 150 — 1450, 
but are shown combined into 3 bins for display purposes. The solid blue line 
is calculated using the mock catalogs and convergence maps. The dashed 
and dotted lines are obtained using the Limber approximation and the red- 
shift distribution from the mock catalogs for the non-linear and linear power 
spectrum, respectively. The upper and lower curves for each model corre- 
spond to a redshift-independent galaxy bias of 1.0 and 1.5, respectively. For 
clarity we only plot models for the optical catalogs. Predictions for the IR 
samples are within 15% of the plotted curves. 

We calculate cross-power spectra between the CMB lens- 
ing maps and the maps of fractional galaxy over-density from 
each of the galaxy catalogs. Given the high source density in 
the galaxy catalogs and the cuts designed to enhance unifor- 
mity in source selection, we expect the noise in these corre- 
lations to be dominated by the lensing reconstruction. Noise 
estimates are obtained by cross-correlating the galaxy maps 
with 50 simulated lensing maps (including realistic noise). 
Results are shown in Figure [3 

Significant correlated signa j(~ 4.2 — 5.3a) is detected for 
every catalog, and zero cross-correlation is ruled out robustly 
in all cases. Additionally, we note a conservative magnitude 
cut is used in this Spitzer analysis — extending the catalogs 
one magnitude fainter increases the correlation to > 8.8<r — 
but, as the properties of this deeper catalog are still being 
characterized, we leave detailed study of this sample to future 
work. 

We test for contamination from Galactic cirrus by cross- 
correlating lensing maps and galaxy maps with pre dictions for 
the Galactic dust emission (Finkbeine r et al.|1999| ) and find no 
correlation. We also place a limit on the systematic bias owing 
to tSZ using the simulated tSZ map described in Section 5. 
We mask an identical density of clusters as in the real CMB 
maps and process this map through the quadratic estimator 
used to create the convergence maps. The resulting maps of 
tSZ leakage are cross-correlated with the mock galaxy density 



maps. The tSZ leakage has a small effect, suppressing the 
correlation by ~ 5%. 

A power-law is fit to each measured cross spectrum: Cl = 
A(L /500) n . In all cases a simple power law is a good fit 
and the best-fit amplitudes are significantly above zero. The 
data are consistent with a single power law describing all data 
sets. The constraints on amplitude and power-law index for 
each cross- spectrum are shown in Figure H] and summarized 
in Table 1 . In Figures [3] and [4] we also plot the results from 
the mock catalogs and convergence maps. Additionally we 
plot the theoretical expectation (calculated using the CAMB 
software package, |Lewis et al. 2000 ) for both the non-linear 
and linear matter power spectrum using the galaxy redshift 
distributions from the mock catalogs. We also constrain the 
redshift-independent bias assuming non-linear evolution, by 
fixing the power law index to the non-linear value and evalu- 
ating the ID amplitude probability distribution relative to the 
b = 1 valu e. These results are co nsistent with previous mea- 
surements ( d e la Torre e t al. 2007]). 

The measured power-law slopes are in rough agreement 
with expectations for large-scale structure that is mainly in 
the non-linear regime. The angular correlation function of 
optical galaxies is reasonably well-fit by a simple power law 
with w(6) oc 0~ - 77 (e.g., Peebles 1975), which corresponds 
to n = —1.23. The fluctuations in the cosmic IR ba ckground 
are also well-fit by a power law with n = —1.2 ([Addison 



2 Signal-to-noise is defined as ^/Ax 2 (0). Ax 2 (0) is the difference in x 2 
from best-fit to a model with zero cross-correlation 



et al.|2012t |Planck Collaboration et al.|20111 [Reich ardt et al. 
201 ip . It is possible that at low L the cross-correlation is prob- 
ing scales that are still in the linear regime, while at high-L 
the clustering power spectrum which dominates has a steeper 
slope, resulting in a good approximation of a power law. The 
data are not yet sufficiently constraining to test this hypothe- 
sis. 

Modeling the signal further is complicated by the lack of 
precise redshift information for each catalog, residual con- 
tamination of the galaxy catalogs by Galactic objects, and in- 
completeness due to source confusion. These are not funda- 
mental limitations, and work is ongoing to both improve and 
better characterize these galaxy catalogs. Even in the pres- 
ence of these limitations, the cross- spectra generally appear 
to agree with the expectations based on both the mock cata- 
logs and theory. 

7. CONCLUSIONS 

We have demonstrated that CMB lensing convergence maps 
are strongly correlated with other tracers of large scale struc- 
ture. These results are derived from a small fraction (< 10%) 
of the now completed 2500 deg 2 SPT survey, and work is 
underway to create convergence maps of the full survey. Re- 
alizing the full potential of the cross -correlation analyses will 
require deep galaxy surveys over the complete survey area. 
The positive correlation with WISE is promising, as all-sky 
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Fig. 4. — Power-law fits to SPT/galaxy density cross-spectra. The 
marginalized 1-d probability distributions for the power-law indices are 
shown projected onto the vertical axis. Also shown are the predictions from 
mock catalogs for optically selected galaxies (star), and predictions using the 
mock redshift distribution and using non-linear (triangle) and linear (square) 
power spectra, assuming a bias of 1 . 

WISE catalogs have recently become available. The SPT- 
Spitzer Deep Field will ultimately cover three times the area 
presented here. This same field is scheduled for observation 
with Herschel-SPIR E and will soon be imaged with SPTpol 
( |Bleem et al.|2012| ), resulting in even deeper lensing maps. In 
the next few years, the Dark Energy Survey and VISTA Hemi- 
sphere Survey are expected to fully cover the SPT survey at 
optical and near-IR wavelengths. 



With existing and planned high signal-to-noise, large-area 
CMB lensing maps, the near future promises rapid evolution 
in the state of CMB lensing analyses. Combining the lens- 
ing maps with multi- wavelength information from upcoming 
galaxy surveys will allow exploration of galaxy formation 
models out to high redshift and of structure growth through 
measurement of the amplitude and evolution of cosmic den- 
sity fluctuations. Careful modeling of the survey galaxy pop- 
ulations, sources of bias in the lensing reconstructions, and 
theoretical expectations for the signal will allow these cross- 
correlation analyses to achieve their full potential. 
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